The Innermost Stable Circular Orbit of Binary Black Holes 
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We introduce a new method to construct solutions to the constraint equations of general relativity 
describing binary black holes in quasicircular orbit. Black hole pairs with arbitrary momenta can be 
constructed with a simple method recently suggested by Brandt and Briigmann, and quasicircular 
orbits can then be found by locating a minimum in the binding energy along sequences of constant 
horizon area. This approach produces binary black holes in a "three-sheeted" manifold structure, as 
opposed to the "two-sheeted" structure in the conformal- imaging approach adopted earlier by Cook. 
We focus on locating the innermost stable circular orbit and compare with earlier calculations. Our 
results confirm those of Cook and imply that the underlying manifold structure has a very small 
effect on the location of the innermost stable circular orbit. 
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I. INTRODUCTION 

Binary black holes are among the most promising 
sources of gravitational radiation for the new genera- 
tion of gravitational wave detectors such as the Laser 
Interferometric Gravitational Wave Observatory (LIGO), 
VIRGO, GEO and TAMA. This has motivated an intense 
theoretical effort to predict the gravitational waveform 
emitted during the inspiral and coalescence of two black 
holes 0. 

Because of the circularizing effects of gravitational ra- 
diation damping, we expect the orbits of close binary 
systems to have small eccentricities. The inspiral of a bi- 
nary black hole system then proceeds adiabatically along 
a sequence of quasicircular orbits up to the innermost 
stable circular orbit (hereafter ISCO), where the evolu- 
tion is expected to change into a rapid plunge and coa- 
lescence 0. The ISCO therefore leaves a characteristic 
signature in the gravitational wave signal, and knowledge 
of its location and frequency is thus very important for 
the prospect of future observations. 

While various approximations may be adequate to 
model the adiabatic inspiral up to the ISCO, it is gen- 
erally expected that only numerical simulations in full 
general relativity can accurately model the dynamical 
plunge and merger and predict the gravitational signal 
from that phase. It is therefore desirable to construct ini- 
tial data for numerical evolution calculations describing 
binary black hole pairs at the ISCO, which adds another 
motivation for determining the location of the ISCO. 

Various approaches have been adopted to locate the 
ISCO in compact binaries, including first order post- 
Newtonian approximations |^] , variational principles , 
second order post-Newtonian methods combined with a 
"hybrid" approach Q], a Fade approximation Q and 
an effective-one-body approach [^,^|, and numerical so- 
lutions to the constraint equations of general relativ- 
ity ||l^,|ll[. Unfortunately, however, the results differ sig- 
nificantly and yet have to show any sign of convergence 
(see Table II below). It would clearly be desirable to un- 



derstand the origin of these differences. In this paper, 
we revisit binary black hole solutions to the constraint 
equations, and evaluate how some of the choices which 
have to be made in this approach affect the location of 
the ISCO. 

Before the constraint equations of general relativity 
can be solved, a background geometry and topology 
have to be chosen. In the conformal-imaging approach 
adopted by Cook [0, a conformally flat (spatial) back- 
ground metric is chosen toget her with a two-sheeted man- 
ifold structure (see Sec. II A ). It has been suggested that 
these choices may affect the location of the ISCO, and 
may explain the difference between these and the more 
recent post-Newtonian results. 

In this paper, we combine the methods of Cook |Tc| ] 
and Brandt and Briigmann fl^ to introduce a new ap- 
proach to constructing binary black holes in quasicircu- 
lar orbit. We follow Cook (lO| and choose a conformally 
flat background metric, but do not assume an inversion- 
symmetry as is done in the conformal-imaging approach. 
This considerably simpl ifies the solution of the momen- 
tum constraint (see Sec. 11 B), and produces binary black 
holes in a three-sheeted manifold structure as opposed to 
the two-sheeted structure in the conformal-imaging ap- 
proach. Moreover, adopting the "puncture" approach of 
Brandt and Briigmann Jl^, the Hamiltonian constraint 
can be solved very easily numerically on R'^ without hav- 
ing to impose boundary conditions on interior bound- 
aries (see Sec. II C ). We locate the ISCO, and find 
that its physical parameters agree very well with those 
found with the conformal-imaging approach of Cook [0 . 
We therefore conclude that the choice of the underlying 
manifold structure has a very small effect on the location 
of the ISCO. Our new approach, which is significantly 
simpler than the conformal-imaging approach, may also 
provide a framework in which the conformal-fiatness as- 
sumption may be relaxed, and its effect on the ISCO be 
evaluated. 

The paper is organized as follows. In Sec. II, we intro- 
duce the basic equations and explain how binary black 
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holes in quasicircular orbit can be constructed. We dis- 
cuss our numerical implementation in Sec. III. In Sec. IV 
we present our results and compare with those from other 
approaches. We briefly summarize in Sec. V. 

II. SETUP OF THE PROBLEM 
A. The Initial Value Problem 

A framework for constructing initial data describing 
binary black holes has been provided by Arnowitt, Deser 
and Misner's 3 + 1 decomposition of Einstein's equa- 
tions ||l^ and York's conformal decomposition JT^ . p^ . 

The 3-1-1 decomposition splits Einstein's equations into 
evolution and constraint equations for the metric 7^ of a 
spatial hypcrsurface S, and the extrinsic curvature ify , 
which describes the embedding of the hypcrsurface S in 
the full spacetime. The physical metric ^ij can now be 
decomposed into a conformal factor ijj and a conformal 
background metric 7^ , 

7»i=V'*7»i- (1) 

It is also convenient to decompose the extrinsic curvature 
Kij into its trace K and a trace-free conformal back- 
ground extrinsic curvature Aij according to 

= i^-^A,, + ^-f,jK. (2) 

The Hamiltonian constraint then reduces to an equation 
for the conformal factor i/j, 

SV^^-i^R-^ij^K^ + ij-^AijA'^ =0, (3) 

and the momentum constraint can be written 

DjA'^ - -i!^fW.jK = 0. (4) 
3 

Here Di is the covariant derivative compatible with the 
conformal background metric, the Laplacian, and R 
is the Ricci scalar. 

Binary black hole initial data cannot be constructed 
uniquely, because the constraint equations of general rel- 
ativity determine neither the background geometry nor 
the topology of the spacetime, both of which have to 
be chosen before the constraint equations can be solved. 
Loosely speaking, these ambiguities correspond to differ- 
ent amounts of gravitational radiation in the initial data 
sets. In this paper, we will follow Cook et.al. 0,^,0 
and choose a flat background geometry, but we will 
choose a three-sheeted topology as opposed to the two- 
sheeted topology of Cook. 

Choosing the background geometry amounts to choos- 
ing the conformal background metric 7^ . Following Cook 
et. al. |l^,|l^,|l^ , we choose the conformal background ge- 
ometry to be flat so that 7ij — fij, where fij is the 



flat metric in a so far arbitrary coordinate system. The 
covariant derivative then becomes the flat-space co- 
variant derivative, and the Ricci scalar R vanishes. We 
will later specialize to cartesian coordinates, 7^^ — Sij, 
for which Di reduces to a partial derivative. We also 
take the hypersurface E to be maximally embedded in 
the spacetime so that K = 0. With these choices, the 
constraint equations simplify to 

VV = -^V'-'iyi'^' (5) 

and 

DjA'^ = 0. (6) 

Note that maximal slicing A' = automatically decou- 
ples the momentum constraint from the Hamiltonian con- 
straint. 

Choosing the topology of the spacetime is less straight- 
forward (compare the discussion in Since we are 
interested in isolated black-holes systems, it is natural 
to assume the hypersurface E to be asymptotically flat. 
Constructing black hole data in vacuum, however, neces- 
sarily involves non-trivial topologies. This can be illus- 
trated by a t — const slice of the Schwarzschild geome- 
try in isotropic coordinates, where every point inside the 
black hole's throat can be mapped into a point outside 
the throat and vice versa. Moreover, such a mapping 
can be accomplished with an isometry, which maps the 
metric into itself, implying that the physical fields at a 
point inside the throat are identical to those at a point 
outside the throat. In particular, the geometry near the 
center is identical to the geometry near infinity. We can 
therefore think of this solution as describing two iden- 
tical, asymptotically flat "universes" or "sheets", which 
are connected by a throat or Einstein- Rosen bridge . 

There is no unique generalization of this topology to 
the case of multiple black holes For two black holes, 
the two throats could either connect to the same asymp- 
totically flat sheet, or else to two separate asymptotically 
flat sheets. The former approach results in a two-sheeted 
topology, the latter in a three-sheeted topology. 

Cook et.al. [ p^Hl^Jl7[| implemented a "conformal-im- 
aging" formalism, which adopts a two-sheeted topology 
together with the additional demand that the two sheets 
are related by an isometry so that their physical fields are 
identical (c/. It has been argued that this choice 

is the "most faithful generalization of the Schwarzschild 
geometry to the case of multiple holes" Moreover, 
the isometry conditions on the throats can be used as 
boundary conditions in numerical implementations, so 
that singularities inside the throats can be eliminated 
from the numerical grid. The computational disadvan- 
tage of this method is that boundary conditions have to 
be imposed on fairly complicated surfaces. In finite dif- 
ference algorithms, this can be accomplished either with 
bispherical or Cadez coordinates ||2^ , designed such that 
a constant coordinate surface coincides with the throat, 
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or else with fairly complicated algorithms in cartesian 
coordinates. Both approaches, together with a spectral 
method, have been compared in Q. 

In this paper, we choose instead a three-sheeted topol- 
ogy and do not assume an inversion-symmetry across the 
throats, which simplifies the problem in two respects. 
The analytical solution to the momentum constraint be- 
comes very simple, since we no longer need to construct 
inversion-symmetric solutions (see Sec. [IB). Moreover, 
the singularities inside the black holes can be removed an- 
alytically using a "puncture" method r ecent ly suggested 
by Brandt and Briigmann |l2| (see Sec. II C). The prob- 
lem can then be solved quite easily on R"^ in cartesian 
coordinates, without having to impose interior bound- 
ary conditions. The only added complication is that one 
now has to locate apparent horizons in the numerically 
constructed hypersurface. 



(10) 



From (0) and (||) we find the that total momentum of 
this solution is P = Pi -I-P2, and the angular momentum 
about the origin of the coordinate system 



J = CiXPi + C2XP2. 



Note that constructing inversion-symmetric solutions 
for multiple black holes in the conformal-imaging ap- 
proach is fairly complicated. There, the components of 
the extrinsic curvature are expressed in terms of an infi- 
nite series of recursively defined quantities (see [pT| and 
Appendix A of px[|). Relaxing the inversion symmetry, 
so that the extrinsic curvature can be written as a sim- 
ple superposition of two solutions (||), therefore greatly 
simplifies the problem. 



B. Solving the Momentum Constraint for Binary 
Black Holes 



For maximally sliced hypersurfaces, the momentum 
constraint decouples from the conformal factor, and ana- 
lytical solutions to (||) can be given. Moreover, for confor- 
mally and asymptotically flat data, the total (physical) 
linear momentum [E3| 



Stt J 00 

and the total (physical) angular momentum 



(7) 



(8) 



can be determined from A''^ without having to solve the 
Hamiltonian constraint (^) (see |Q). 

Analytical solutions to the momentum constraint (|^) 
describing single boosted or spinning black holes have 
been given by Bowen and York p^-|2^. A solution A^^ , 
describing a single black hole at the coordinate location 
C with linear momentum P is given by 



^' c ^ 



(9) 



Here re = — C*|| is the coordinate distance to the cen- 
ter of the black hole and n\~, — (x' — C'^)/r is the normal 
vector pointing away from that center. Additional terms 
have to be added in the conformal-imagine approach for 
an isometry condition to hold across the throat. Note 
that we have only included linear momentum terms in 
this expression, and that we are therefore restricting our 
analysis to non-spinning black holes. 

Since the momentum constraint (|^) is linear, we can 
construct binary black hole solutions by superposition of 
single solutions 



C. Solving the Hamiltonian Constraint for Binary 
Black Holes 

Solutions to the Hamiltonian constraint (|^) can be 
constructed by generalizing the Schwarzschild solution 
in isotropic coordinates for a static {i.e. Aij = 0) and 
spherically symmetric black hole at coordinate location 
C, 



(12) 



(note that asymptotic flatness demands t/i ^ 1 as r — > 
00). Solutions describing multiple static black holes can 
be constructed by adding contributions M./{2rc) for 
each black hole. To establish an inversion-symmetry, ad- 
ditional terms would again have to be added pof . 

In the "puncture" method suggested by Brandt and 
Briigmann , a general nonstatic solution to the Hamil- 
tonian constraint is written as a sum of the static, an- 
alytic contribution plus a term correcting for finite Aij. 
Adopting their notation, we write 



tp = u + 



1 



where a is defined by 
1 



Ml 



a 2r, 



M2 
2rc2 



(13) 



(14) 



The Hamiltonian constraint then becomes an equation 
for the correction term u 



V^u = -/3(l + au)"^ 
where we have abbreviated 



(15) 



(16) 



For asymptotic fiatness, we impose a Robin boundary 
condition d{r {u —l))/dr = at large distances from the 
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black holes. The existence and uniqueness of solutions 
u on R'^ has been established in [|l2|. The beauty of 
this approach is that the poles at the center of the black 
holes have been absorbed into the analytical terms. The 
corrections u are regular everywhere and can be solved for 
very easily on a simple computational domain, without 
having to impose boundary conditions on the throats. 

Once the conformal factor ^ has been determined, the 
ADM mass of the solution can be found from 



E 



1 
1 

2^ 



1 



a I 2tt 



= TWi + M2 + 7?- / /3 (1 + auy^dV. 

l-K J 

Note that the integral extends over all space. 



(17) 



D. Constructing Equal- mass Binary Black Holes in 
Quasicircular Orbit 



two apparent horizons, which is a very good approxima- 
tion to the shortest proper separation between the two 
horizons. 

Quasicircular orbits can then be found quite easily 
(see [^) by computing the effective potential E(, as a 
function of separation I along a sequence of constant 
black hole mass M and angular momentum J and lo- 
cating turning points 



dEt 
dl 



0. 



(22) 



A minimum corresponds to a stable quasicircular orbit, 
while a maximum corresponds to an unstable orbit. For 
a quasicircular orbit, the binary's orbital angular velocity 
f2 as measured at infinity can then be determined from 



dEt, 
dJ 



(23) 



MI 



(see footnote for a Newtonian illustration). 



We now specialize to equal mass black holes with M = 
Ml = A42- In the center-of-mass frame of the binary 
system, we have 



Pi 



(18) 



Binaries in quasicircular orbit should furthermore satisfy 
P • C = 0, where we have defined 



C = Ci — C2 



(19) 



Without loss of generality, we can then take C to be 
aligned with the 2-axis, P to be aligned with the x-axis, 
and place the origin of the cartesian coordinate system 
at the center between the two black holes. 

The problem has now been reduced to a three-dimen- 
sional parameter space with the free parameters A4, C = 
||C|| and P = ||P||. For every configuration, we compute 
several physic al q uantities. We determine the total ADM 
mass E from (O) and the total angular momentum J = 
Jy = PC from (|11|) . Since we have restricted our analysis 
to non-spinning black holes, the mass of each individual 
black hole can be identified with the irreducible mass 



M = Mir 



A 

16^ 



1/2 



(20) 



where A the proper area of the black hole's apparent 
horizon |^ . We now define the effective potential as the 
binding energy 



Eh^E- 2M. 



(21) 



Lastly, we compute the proper separation / between the 
two horizons along the line connecting the centers of the 



III. NUMERICAL IMPLEMENTATION 

We adopt a finite difference approach to solve eq. ( p^ ) 
in cartesian coordinates. The numerical code is imple- 
mented in a parallel, distributed memory environment 
using DAGH software ||2^], and the Laplace operator 
in (15) is inverted using PETSc software |Q. We lin- 
earize eq. (|l^) and iteratively solve for corrections to 
approximate solutions until convergence to a desired ac- 
curacy has been achieved. Since the components of the 
extrinsic curvature ( p^ are either symmetric or antisym- 
metric across the coordinate planes x = 0, y — and 
2: = 0, it is sufficient to solve the Hamiltonian constraint 
in only one octant. 

In addition to verifying second order convergence of 
our code, we have performed tests in the linear regime 
by comparing with the linear analytic solution for black 
holes boosted towards each other (eq. (16) in ||l2|), and 
in the nonlinear regime by comparing with the "A2B8" 
dataset jl^], for which values of the ADM mass in a 
three-sheeted manifold structure have been given in Ta- 
ble I of |l^. Note, however, that those masses have er- 
roneously been calculated for black hole spins with signs 
opposite to those given in that table, "Si, 2 = —81,2" [ pl| . 

Given a solution for a set of parameters Ai, C and 
P, we can locate an apparent horizon and determine the 
black hole mass AI from the horizon's proper area using 
the algorithm described in . This algorithm expresses 
the location of a closed surface in terms of symmetric 
trace-free tensors, and varies the expansion coefficients 
until an outer-most trapped surface has been found. We 
found that for the horizons in this problem, which are 
fairly spherical and only very mildly deformed, an ex- 
pansion up to order £ = A is adequate. 
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Before constructing sequences of constant black hole 
mass M , it is convenient to rescale all variables with re- 
spect to that desired value of the black hole mass M . We 
introduce the sum of the black hole masses 



Ml + M2 = 2M 



and the reduced mass 



M1M2 
Ml + M2 



M 

Y' 



(24) 



(25) 



and define the dimensionless parameters M = M./m 
C = C/m and P = We also rescale the angu- 

lar momentum and angular velocity as J = J/ and 
f2 = mQ, and identify the dimensionless effective poten- 
tial with the rescaled binding energy Eh = Eh/ii. 

Sequences of constant J can now be constructed by set- 
ting P = J/C for a set of different values of C, and by it- 
erating over M until the (dimensionless) numerical value 
of the black hole mass M = Mnum/M has converged to 
unity within a desired accuracy for each C |]33|] . 

For typical cases of interest, the difference between the 
ADM mass E and the sum of the black hole masses m 
is quite small. According to (pT|), the binding energy is 
therefore much smaller than those masses, and its relative 
numerical error much larger. In order to reliably locate a 
minimum in the binding energy, we therefore have to de- 
termine the masses to very high accuracy (c/. |p^). More- 
over, numerous models have to be calculated to construct 
a sufficient number of sequences over a sufficient range of 
separations, which is only affordable for a fairly moder- 
ate maximum grid size. For a uniform, cartesian grid 
of a given size, a compromise then has to be found be- 
tween extending the computational grid to large enough 
distances and sufficiently resolving the individual black 
holes. 

The location of the outer boundary of the computa- 
tional grid affects the results through the Robin bound- 
ary condition on u, which is correct only asymptotically, 
and the energy integral ([T7|), which should extend over 
all space. The latter effect can be improved by extend- 
ing the energy integral beyond the numerical grid, where 
a and (3 can be evaluated analytically, and where u can 
be estimated from its value on the surface of the com- 
putational grid and its 1/r falloff. In the top panel of 
Fig. 1, we show the binding energy for different locations 
of the outer boundary Zo^t = ^out/"^ for a typical con- 
figuration of interest (C ~ 2.5, J — 3.0) with a fixed grid 
resolution h = h/m = 0.03125. For all calculations pre- 
sented in this paper we use Xout = ^out = ^out/2. For 
the dashed line, only contributions to the ADM mass 
from inside the computational grid have been taken into 
account, and for the solid line we have expanded the vol- 
ume for the energy integral by a factor of six in each 
dimension. Obviously, the latter converges much more 
rapidly and yields a more accurate value for all locations 
of the outer boundary. 
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-OA r 
-0,12 r 
-0.14 




0.02 0.03 0.04 0.05 0.06 



FIG. 1. The binding energy Et for C = 2.5 and J = 3.0 for 
different locations of the outer boundary Zout (at constant 
grid resolution h = 0.03125, top panel) and different grid 
resolutions h (at constant outer boundary Zout = 6, bottom 
panel). The dashed line only includes contributions to the 
ADM mass from inside the computational grid, and the solid 
line denotes the corrected value (see text). 



The resolution of the individual black holes affects the 
accuracy with which their apparent horizons can be lo- 
cated, and hence the accuracy of their masses M. The 
effect of the resolution on the binding energy is demon- 
strated in the bottom panel of Fig. f , where we show Et, 
for different grid resolutions for the same configuration, 
this time with the outer boundary fixed at Zout — 6. 

From Fig. 1, we find that the binding energy Eb can 
be determined to within at most a few percent error for 
Zout = 6 and h = 0.03125, corresponding to a numerical 
grid of size 96 x 96 x 192. The iteration to construct 
one model then takes approximately 3 CPU hours on the 
NCSA Origin2000, which makes an extensive survey of 
parameter space affordable. We use these grid specifica- 
tions for all results presented in the following section. 



IV. RESULTS 

In Fig. 2, we show contour plots of the conformal factor 
-0 for a configuration close to the ISCO {C = 2.25 and 
J ~ 2.95). The apparent horizons, marked by the thick 
dashed lines, are dragged along by the black holes and 
lag slightly behind in their counter-clockwise orbit. This 
effect has been discussed for single boosted black holes 
in Q . Note that we compute the proper separation I = 
l/ra between the horizons along the line connecting the 
centers of the apparent horizons. This is a coordinate- 
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FIG. 2. Contours of the conformal factor i/> for a configu- 
ration close to the innermost stable circular orbit (C — 2.25 
and J = 2.95). The contours (solid lines) logarithmically span 
the interval = 1 and tp = 9.2. Note that the apparent hori- 
zons, marked by the thick dashed lines, are not concentric 
with the contours of the conformal factor. Instead, they are 
dragged along by the black holes and lag slightly behind in 
their (counter-clockwise) orbit. 



I 


Et 


J 


n 


P/a 


C/a 


4.8 


-0.092 


2.95 


0.18 


1.7 


5.9 


4.880 


-0.09030 


2.976 


0.172 


1.685 


5.91 



TABLE I. Comparison of our results for the innermost sta- 
ble circular orbit (top line) with those of Cook (bottom 
line). 



dependent quantity, but a very good approximation to 
the (coordinate-independent) shortest proper separation 
between the horizons. 

We now construct sequences of constant angular mo- 
mentum for various values of J, and plot the effective 
potential Eb along these sequences as a function of I 
in Fig. 3. A minimum in the effective potential cor- 
responds to a stable quasicircular orbit. The bold line 
connecting these minima in Fig. 3 represents a sequence 
of quasicircular orbits. This sequence terminates at the 
ISCO, where the adiabatic, quasicircular inspiral of the 
two black holes is expected to change into a rapid plunge 
and merger j2j . Since this transition leaves a characteris- 
tic signature in the gravitational wave signal, the knowl- 
edge of the location and frequency of the ISCO are of 
great importance for future observations with the new 
generation of gravitational wave detectors. 

In Table I, we list our results for the physical parame- 



4 5 6 7 8 
I 

FIG. 3. The effective potential Eb as a function of proper 
separation / for the following values of the angular momen- 
tum J: 2.9, 2.92, 2.94, 2.95, 2.96, 2.97, 2.98, 3.00, 3.02, 3.06, 
3.10, 3.15, 3.20, 3.25 (from bottom to top). Quasicircular or- 
bits correspond to minima in the effective potential. The bold 
line connects these minima and represents a sequence of qua- 
sicircular orbits. This sequence terminates at the innermost 
stable circular orbit. 

ters of the ISCO (top line) and compare with the results 
of Cook (bottom line). We tabulate the proper sep- 
aration between the horizons I, the binding energy Eh, 
the angular momentum J, and the angular velocity n 
as well as the the linear momentum of each black hole 
P/a and the coordinate separation C/a, where a is the 
(coordinate) radius of the black holes. The latter is not 
well-defined in our calculation, but since the black holes 
are nearly spherical it is very reasonable to estimate an 
average radius from the £ = monopole term in the mul- 
tipole expansion for the horizon. 

We conclude that all quantities agree fairly well with 
those of Cook ||l0|] within our estimated numerical error of 
a few percent. As the most significant deviation, we find 
that in our calculation the binary is slightly more tightly 
bound at the ISCO, and correspondingly has a slightly 
larger angular velocity. However, even these quantities 
differ by less than 5%, which may be caused by nu- 
merical effects. We conclude that the choice of the un- 
derlying manifold structure has a very small effect on the 
location of the ISCO. 

In Table II we compare the binding energy, the angular 
momentum and the angular velocity at the ISCO from 
various different calculations. For a test particle orbit- 
ing a Schwarzschild black hole, the ISC O can be located 
analytically, which yields Eb = \/S/Q ~ 1 ^ -0.0572, 
J = 2^3 - 3.464, and n = 1/6^/2 - 0.0680. Clark 
and Eardley l5l adopted a first order post-Newtonian ar- 
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Reference 



Eh 



J 



Schwarzschild 
CE |l 
BD g 
KWWjl 
Cook |o| 
BCSST|l| 
DIS i) 
BD gl 
This work 



-0.0572 

-0.1 

-0.65 

-0.0378 

-0.09030 

-0.048 

-0.0653 

-0.06005 

-0.092 



3.464 

3.3 

0.85 

3.83 

2.976 

3.9 

3.40 
2.95 



0.068 
2 

0.0605 

0.172 

0.06 

0.0885 

0.0734 

0.18 



TABLE II. Comparison of the binding energy Et, the an- 
gular momentum J and the angular velocity H for the inner- 
most stable circular orbit from various different calculations: 
a test particle orbiting a Schwarzschild black hole, Clark and 
Eardley (CE) ^, Blackburn and Detweiler (BD) [||, Kidder, 
Will and Wiseman (KWW) §, Cook 0, Baumgarte et.al. 
(BCSST) 0, Damour, Iyer and Sathyaprakash (DIS) g, 
Buonanno and Damour (BD) and the results from this 
paper. The results of Baumgarte et.al. |llj are for an n = 1 
polytrope binary of compaction M/R = 0.2. Naively extrap- 
olating to M/R — 0.5 yields values very close to our results 
for binary black holes. 



gument to approximately estimate the location of the 
ISCO. We list their values for nonrotating neutron stars. 
Blackburn and Detweiler j^] adopted a variational prin- 
ciple and assumed a periodic solution to Einstein's equa- 
tions. At the ISCO, for which the approximations of this 
approach fail, they find extremely tightly bound bina- 
ries. Kidder, Will and Wiseman adopted a second or- 
der post-Newtonian approximation together with a "hy- 
brid" approach and found an extremely weakly bound 
ISCO. However, several authors have cast doubt on the 
robustness and consistency of the hybrid approach |^,|| . 
Baumgarte et.al. [|ll| constructed fully relativistic mod- 
els of corotating binary neutron star in quasiequilibrium, 
albeit assuming conformal flatness, and found that the 
ISCO depends on the compaction of the neutron stars. 
In Table I we list their results for n = 1 polytropes of 
compaction M/R = 0.2, where M and R are the mass- 
energy and areal radius which the stars would have in iso- 
lation. Naively extrapolating their results to M/R = 0.5 
yields values which are very similar to our result for bi- 
nary black holes. Damour, Iyer and Sathyaprakash H] 
combined a second order post-Newtonian approximation 
with a Pade approximation, which yields a slightly more 
tightly bound ISCO than that for a test-particle orbit- 
ing a Schwarzschild black hole. A similar result is found 
by Buonanno and Damour, who combine a second or- 
der post-Newtonian approximation with an effective-one- 
body method MM. 



V. SUMMARY AND DISCUSSION 

Since an accurate knowledge of the ISCO is very impor- 
tant for possible future gravitational wave observations, 
it is very unsettling that different approaches to com- 
puting the ISCO lead to very different results (compare 
Table II). One of these approaches, namely constructing 
binary black hole solutions to the constraint equations of 
general relativity, involves choosing the background ge- 
ometry and topology, and it would be very desirable to 
know how much the results depend on these choices. 

In this paper, we introduce a new method to construct 
solutions to the constraint equations of general relativ- 
ity describing binary black holes in quasicircular orbit. 
We combine the approaches of Cook and Brandt 
and Briigmann to construct binary black holes in 
a three-sheeted manifold structure, as opposed to the 
two-sheeted topology in the conformal-imaging approach 
adopted by Cook |10|. We locate the ISCO and find that 
its physical parameters are very similar to those found by 
Cook Our results confirm those earlier results and 
imply that the underlying manifold structure only has 
a very small effect on the ISCO. The latter is perhaps 
not entirely surprising, since it reflects the fact that the 
strength of the imaged poles in the conformal-imaging 
approach is smaller than the strength of the poles them- 
selves [|o|. 

Our new approach is considerably simpler than the 
conformal-imaging approach of Cook [Q. The analytic 
solution to the momentum constraint simplifies because 
no inversion-symmetric solutions have to be constructed, 
and the numerical solution to the Hamiltonian constraint 
simplifies because we can adopt the "puncture" method 
of Brandt and Briigmann jl2|. In particular, we can solve 
the Hamiltonian constraint in cartesian coordinates on 
R^ without having to impose interior boundary condi- 
tions. One disadvantage of our approach is that the ap- 
parent horizons have to be located numerically, which we 
do with the algorithm developed i n p2[ . 

In this paper, we follow Cook [|10[ and choose a con- 
formally flat background metric. Accordingly, we cannot 
address the dependence of the ISCO on the choice of the 
background geometry. However, since our new method 
is significantly simpler than the conformal-imaging ap- 
proach, it may provide a useful framework to relax the 
assumption of conformal flatness and to construct binary 
black holes in quasicircular orbit for more general back- 
ground geometries. 
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